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Abstract. In order to put strict observational constraints on the evolutionary status of the magnetic chemically 
peculiar stars (CP2) of the upper main sequence, we have investigated a well established sample of galactic field 
CP2 objects within a radius of 200 pc from the Sun in the (X,Y) plane. In total, 182 stars with accurate parallax 
measurements from the Hipparcos satellite were divided into Si, SiCr and SrCrEu subgroups based on classification 
resolution data from the literature. Primarily, it was investigated if the CP2 phenomenon occurs at very early 
stages of the stellar evolution, significantly before these stars reach 30 % of their life-time on the main sequence. 
This result is especially important for theories dealing with stellar dynamos, angular momentum loss during the 
pre- as well as main sequence and stellar evolutionary codes for CP2 stars. For the calibration of the chosen 
sample, the well-developed framework of the Geneva 7-color and Stromgren uvbyfi photometric system was used. 
We are able to show that the CP2 phenomenon occurs continuously at the zero age main sequence for masses 
between 1.5 and 4.5 Mq. The magnetic field strengths do not vary significantly during the evolution towards the 
terminal age main sequence. Only the effective temperature and magnetic field strength seem to determine the 
kind of peculiarity for those stars. We found several effects during the evolution of CP2 stars at the main sequence, 
i.e. there are two "critical" temperatures where severe changes take place. There is a transition between Si, SiCr 
and SrCrEu stars at 10000 K whereas a significant decrease to almost zero of evolved SrCrEu objects with masses 
below 2.25 Mq at 8000 K occurs. These conclusions have to be incorporated into models that simulate the stellar 
formation and evolution of stars between 1.5 and 4.5 Mq in the presence of strong magnetic fields. 
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^ ■ 1. Introduction — The dynamo theory based on the existence of a con- 

i , , . temporaneous dynamo operating in the convective core 

In Pohnl et al. (2003) we nave shown that young magnetic r ,. m , . , , , 

, . \. ' i i- 7. i °i the magnetic stars, lhis means that the magnetic 

2 . chemically peculiar stars (CP2 neraiter, according to the n , , , , , , , , . , 

, r .„ . ° holds should only appear at evolutionary stages De- 
notation by Preston 1974 oi the upper mam sequence ex- , n i, • i ■ /rrAnio 

' , . twccn the zero and terminal age mam sequence (/AMa 

ist m young open clusters (ages between 10 and 140 Myr) and TAMg hereafter) preferabl closer to the TAMS. 

of our Milky Way. Our result received significant support , ., ,, , , . „ ,, 

J . J ,_ , , _„-„ s . , . _ — ihe fossil theory has two variants: the magnetic held 

by the detection Uagnulo ct al. 2003 of a 14.5 kO mag- .,, n , , . ,. „ „ . . 

, . n i , r im ° „ • ^t^,^, ^r-J. , • , , is either the slowly decaying relic ot the rrozen-m m- 

netic held tor HO 66318 m NOC 2516 which was also an , „ ,. „ , , „ , ,. . ., 

„ , , , „„, , , . „„. tcrstellar magnetic field or ol trie dynamo acting m the 

object of the study by Pohnl et al. 2003). 5 , „ J ° ,,. 

. n , , „ , , . pre-mam sequence phase. However, according to this 

Ihe strong stellar magnetic holds oi these omccts to- , , ,. „ , , , , , n , 

..,,.„. ° . . .. model, strong magnetic helds should be present when 

gether with the diffusion oi chemical elements depending , 74 , , c 

. . a star readies trie ZjAIVIo. 
on the balance between gravitational pull and uplift by the 

radiation field through absorption in spectral lines mani- c nj . , . , , ,. £ , , , ,. , 

off several divergent interpretations ot tire evolutionary sta- 

fests in strong overabundances for heavy elements such as , jj_i.--i-j.-r j.- c uj.i • cr^r,n 

° J tus and their implication for magnetic held theories of CP2 

silicon, chromium, strontium and europium compared to starg can be found ^ ^ literature . We have summarize d 



the Sun 

The most important aspect for CP2 phenomenon is 
the origin of the global stellar magnetic fields for those 
objects. Two theories have been developed in this respect 
(Moss 1989): 



them in more detail in Pohnl et al. (2003). 

Our main goal for this paper is to investigate the evo- 
lutionary status for the CP2 group in the solar neighbor- 
hood and to answer the question raised by Hubrig et al. 
(2000) who stated that 1.) The distribution of CP2 stars of 
Correspondence to: ernst.paunzen@univie.ac.at masses below 3Mq in the Hertzsprung-Russell-Diagram 



2 



H. Pohnl, E. Paunzen, H.M. Maitzen: On the formation and evolution of magnetic CP stars 



-50 



* * * 



* 

• • • 



* • 



o 

* o 



* *★ 
• • * * 

★•5 Q* • j- 1 



* 
* 



o 



• * 

* o 



• Si 




o SiCr 




* SrCrEu 




* * 




o * 


• o 






* o * °*c 



.** 



° J5 

* * . * 



-200 -150 -100 -50 50 100 150 200 

X[pc] 

Fig. 1. The spatial distribution of the investigated sam- 
ple of CP2 stars divided into Si (filled circle), SiCr (open 
circle) and SrCrEu (asterisk) objects. 



(HRD hereafter) differs from that of the "normal" stars in 
the same temperature range at a high level of significance: 
magnetic stars are concentrated toward the center of the 
main sequence band and 2.) Magnetic fields appear only in 
stars that have already completed at least approximately 
30% of their main sequence life-time. 

We have investigated three subgroups (Si, SiCr and 
SrCrEu) of CP2 stars which are unambiguously classified 
and have accurate parallax measurements. This sample in- 
cludes 182 objects within a radius of 200 pc from the Sun 
in the (X,Y) plane. The astrophysical parameters of these 
objects were derived taking into account the peculiar char- 
acteristics and all results have been carefully compared 
with already published values. 

2. Program stars 

The target selection for this analysis is based on two main 
criteria: 1.) CP2 stars with known spectral characteris- 
tics and 2.) Accurate parallax measurements available. 
We have taken the lists of Renson (1991), Gomez et al. 
(1998), Renson & Catalano (2001), Bychkov et al. (2003), 
Skiff (2003) and Paunzen et al. (2005) to search for well- 
established CP2 stars. We have only included stars that 
are unambiguously classified within these references. 

From this first list, objects were selected for which the 
absolute error of the parallax is less than 17.5% which is 
the upper limit for applying a correction as suggested by 
Lutz & Kclkcr (1973). Our final list includes 182 objects 
which were divided into Si (61), SiCr (34) and SrCrEu 
(87) stars. These are the three subgroups of the CP2 ob- 
jects that are described in more detail by Preston (1974) 



and Stepien (1994). The phenomenology is due to classifi- 
cation resolution spectra in which different lines of several 
elements are more pronounced than in normal objects of 
the same effective temperature. Even if these subgroups 
are investigated with high resolution spectra, the charac- 
teristic abundance pattern is still different (Wolff 1983). 

Figure ^ shows the spatial distribution of our sample. 
All objects are within a radius of 200 pc from the Sun in 
the (X,Y) plane. 

3. Determination of the evolutionary status for 
the CP2 stars 

In the following subsections we describe the method to 
derive the log L/Lq and logT c ff values for our program 
stars in more detail. Special care was taken to include 
and correct the effects of the spectral peculiarities when 
deriving the age as well as mass of the individual stars. 

3.1. Reddening and metallicity 

Our program stars are located within a radius of 200 pc 
of the Sun in the (X,Y) plane (Fig. fjj for which the in- 
terstellar reddening is not negligible any more. The most 
common way to derive the reddening for single stars is 
to apply standard relations within the Stromgren uvbyfi 
photometric system (Crawford 1975, 1979; Hilditch et al. 
1983). These relations have to be used with caution when 
applied to CP2 stars because all calibrations are primar- 
ily based on the f3 index. Due to variable strong magnetic 
fields, the (3 index can give erratic values (Catalano & 
Leone 1994). Another effect taken into account is a "blue- 
ing" effect which manifests in bluer colors due to stronger 
UV absorption than in normal type stars which can be er- 
ratically interpreted as strong reddening (Adelman 1980). 

We have therefore carefully checked the reddening val- 
ues from the Stromgren uvbyfi photometric system with 
other sources. 

As a first approximation, we have used a "standard" 
value of 0.25 mag kpc -1 derived from open clusters (Dutra 
& Bica 2000). These values were compared with the red- 
dening laws listed by Chen et al. (1998). As the authors 
note, the reddening for very close objects can be overes- 
timated. We found that for several objects the reddening 
derived from Chen et al. (1998) has to be multiplied by a 
factor of 0.7 to be comparable to the other sources. 

As a final step, we have compared and averaged the 
values from all sources. As expected, many stars exhibit 
no reddening at all. 

The metallicity has a severe influence on the determi- 
nation of the age and mass but only marginally affects 
the effective temperature (Glagolevskij 1994). Classical 
CP2 stars are true Population I objects with intrisic so- 
lar abundances. Their peculiarities are restricted to the 
surface only (Nishimura et al. 2004). It is therefore justi- 
fied to use solar isochrones for the determination of the 
age and mass (Gomez et al. 1998), but what metallic- 
ity should be chosen for the solar neighborhood? Schaller 




Fig. 2. The Hertzsprung- Russell-Diagram for our sample of CP2 stars, the symbols are the same as in Fig. ^ The 
evolutionary tracks from 1.5 to 4Mq, lines of equal t/thr, the ZAMS as well as TAMS are for [Z] = 0.016 taken from 
Schaller et al. (1992) and Schaerer et al. (1993). The typical errors are plotted in the upper right corner. 



et al. (1992) give [Z] = 0.0188 whereas Allende Prieto et 
al. (2004) found that the solar neighborhood would be 
"metal- weak" compared to the Sun by 0.06 to 0.11 dex. We 
have therefore used a value of [Z] = 0.016 for the isochrones 
which were interpolated within the grids by Schaller et al. 
(1992, [Z] =0.02) and Schaerer et al. (1993, [Z] =0.008). 

3.2. Effective temperature and luminosity 

For the calibration of the effective temperature, we have 
made use of the Geneva 7-color photometric system 
(Golay 1980). All photometric measurements were taken 
from the General Catalogue of Photometric Data (GCPD, 
Mermilliod et al. 1997). 

The approach by Hauck & Kunzli (1996) was followed. 
For hot CP2 stars (T e // >9500K), the effective tempera- 
ture is given as 

T ef f = -230 + 0.941 • T(X, Y) 

where T(X, Y) has to be calibrated within the grids by 
Kunzli et al. (1997). For cooler objects, we have used the 
calibration by Hauck & North (1993) which is based on 



(B2 — G)q. For the few cool objects without available 
Geneva colors, we have used the Stromgren uvbyft pho- 
tometric calibration by Napiwotzki et al. (1993) which is 
based on the reddening free \c\] index 

5400/T e// = 0.2489 + 0.2698- [ci]. 

A comparison of the results for both calibrations of other 
CP2 stars in the same temperature range shows no system- 
atic trend which is in line with the results by Glagolevskij 
(2002). 

The estimation of the errors for such a calibration has 
been widely discussed in the literature (Smalley & Kupka 
1997). Hauck & Kunzli (1996), for example, list a statisti- 
cal error for the photometric calibration of a large sample 
of normal type stars of ±386 K depending on the tem- 
perature range. North (1998) investigated especially CP2 
stars and found an error of the mean of 4.4% and 3.4% 
(A logTcff = 0.019 and 0.014) for objects hotter and cooler 
than 9500 K, respectively. We were not able to compare 
our calibrated effective temperatures with those used by 
Gomez et al. (1998) because the latter are not available. 
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Fig. 3. A detailed statistical analysis of our sample. We have divided the sample into Si (upper left), SiCr (upper 
right) and SrCrEu (lower left) subgroups and investigated the incidence of these objects depending on their masses. 
The lower right panel shows the comparison of all stars divided into subgroups with masses higher (denoted with a 
pattern) and lower than 3 Mq according to the investigated sample by Hubrig et al. (2000). The Spearman rank order 
correlation coefficient and the % 2 test for these two samples give both a significance level of about 25% for a possible 
correlation. However, these tests are only based on five data points. 



The luminosity calibration for CP2 stars is also af- 
fected by stronger UV absorption which results in correc- 
tion factors Sbc oi u P to 0.3 mag for the estimation of 
bolometric magnitudes for this group (Lanz 1984, North 
1998). The luminosity can be calculated as: 

log L/Lq = 0.4 • (4.72 - M v - BC + S BC ) 

with 

M v = m v + 5 • log 7r + 5 - 3.1 • E(B - V). 

The values for the bolometric correction BC and the stan- 
dard solar value were taken from Flower (1996). Finally, 
we have corrected for the "Lutz-Kclkcr effect" (Lutz & 
Kclkcr 1973) which corrects for the bias in the absolute 
magnitude of a star as estimated from its trigonometric 
parallax. These corrections are rather small compared to 
the error of the parallax measurements from Hipparcos at 
distances up to 200 pc from the Sun. 



3.3. Age and mass 

We have used interpolated isochrones with [Z] = 0.016 
from Schaller et al. (1992) and Schaerer et al. (1993) for 
the determination of the ages and masses for the program 
stars. 

The location of a star within the HRD is unambigu- 
ously defined by the hydrogen concentration in the core 
Xq and the stellar mass. The interpolation of these two 
parameters within the isochrones has been described in 
Pohnl et al. (2003). 

The knowledge of Xq allows us to determine the rel- 
ative age t/thr, where trr denotes the time a star is 
on the main sequence. This is because log L/Lq and 
logM/MQ are, at any stage of the stellar evolution, well 
correlated with Xq Schaller et al. (1992) list value of 
X c = 0.68 and 0.03 for the ZAMS (t/t hr «0) and the 
TAMS (t/thr = 1), respectively. These values are almost 
independent of the stellar mass. 
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Fig. 4. The logi? e versus t/thr as well as \ogT e ff diagrams for our program stars. The magnetic field strengths were 
taken from Bychkov et al. (2003). 



Our quadratic interpolation within the isochrones re- 
sults in: 

t/thr = 0.015 + 2.714 • (0.68 - X c ) - 1.845 • (0.68 - X c f . 

The statistical error for the complete sample, only taking 
into account the errors within the interpolation grid, is less 
than ±5 % for the relative age. The errors for individual 
objects can be larger due to the errors of the effective tem- 
perature (error of the statistical mean is less than 4.4%, 
see Pohnl et al. 2003) and the luminosity (determined be 
the parallax error which is always smaller than 17.5%). 
But throughout this paper, only statistical samples are 
used instead of values for individual CP2 stars which im- 
plies that the statistical error of less than ±5 % for the 
relative age is justified. 

4. Results 

In Fig.[2]the location of the program stars within the HRD 
is shown. We have included evolutionary tracks from 1.5 
to 4M , lines of equal t/t H r, the ZAMS and TAMS for 
[Z] =0.016 and the typical error bars. At first sight it is 
evident that the CP2 stars occupy the whole width of the 
main sequence which means that they are luminosity class 
V objects. 

We performed a detailed statistical analysis of the rel- 
ative ages t/thr for the three samples of Si, SiCr and 
SrCrEu objects. For this purpose we further divided those 
samples according to the masses into four subsamples 
ranging from below 2 Mq to masses higher than 4 Mq . 

Figure OH shows the result graphically. The different 
groups were binned with bin sizes of 0.2 for t/thr. A 



smaller bin size is not justified because of the overall errors 
and the additional introduction of small number statistics. 
From Figs. [21 to 0] we are able to conclude: 

— The CP2 phenomenon occurs already at the ZAMS. 
The percentage of all investigated CP2 stars with 
t/thr < 0.2 is 14.4% as well as 20.8% for objects below 
and above 3Mq which accounts for 16% (or 29 stars) 
of the complete sample. 

— The magnetic field strengths do not vary significantly 
during the main sequence evolution. 

— The age distribution of objects below (no clear 
trend with the age) and above (clear peak at 
0.6 < t/thr < 0.8) 3M is statistically different. The 
Spearman rank order correlation coefficient and the x 2 
test (Rees 1987) for these two samples both give a sig- 
nificance level of about 25% for a possible correlation. 
This result has to be taken with caution because only 
five data points are include in this analysis. 

— The temperature and the magnetic field strength are 
the main parameters that determine the manifestation 
of the chemically peculiarity, i.e silicon - chromium - 
strontium and europium (from hotter to cooler tem- 
peratures) 

— Si stars: they can be found within the whole mass 
range, but more so at young relative ages with higher 
masses. Their incidence seems to decrease with age ex- 
cept for a peak at between 0.6 < t/thr < 0.8. We are 
not able to decide whether this is due to poor number 
statistics or a physical effect. 

— SiCr stars: most of these objects have masses be- 
tween 3 and 4Mq. They seem to be transition ob- 
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jects between Si and SrCrEu stars. There is, for exam- 
ple, an apparent edge at [logT e //,log-B e ] defined by 
[4.0,2.5] and [4.15,3.5] which clearly separates Si and 
SiCr stars. The age distribution clearly peaks between 
0.4<t/t H r<0.6. 
— SrCrEu stars: these are the low mass CP2 stars 
(M < 3Mq) with an upper limit of the effective tem- 
perature between 9500 and 10000 K. There is a clear 
negative correlation between the incidence and the age. 
However, at least 24% of all objects have relative ages 
below 0.4. 

The question arises if the age distribution of CP2 stars 
with masses below 3 Mq is significantly different from that 
of normal type objects in the same spectral range as sug- 
gested by Hubrig et al. (2000). They have used all appar- 
ent single type objects of the Bright Star Catalogue as a 
comparison which resulted in 412 objects. However, their 
CP sample includes only 33 magnetic stars which is most 
certainly biased by lower rotational velocities than for the 
normal type stars (Abt & Morell 1995). The comparison of 
such widely different samples has to be treated with cau- 
tion. We have chosen to follow the approach by Paunzen 
et al. (2002) who used the sample of the program stars to 
choose a random sample of normal type objects by match- 
ing the effective temperature distribution. They found for 
stellar masses below 2.7 Mq a peak of the relative age 
distribution for 0.6 < t/thr < 0.9 (see Figure 2, Paunzen 
et al. 2002). We found exactly the same behaviour for a 
random sample of normal type objects within a radius 
of 200 pc of the Sun in the (X,Y) plane which matches 
the sample of CP2 stars with masses below 3.0 Mq. We 
therefore conclude that the different age distribution of 
low mass CP2 stars as discussed by Hubrig et al. (2000) is 
due to a bias of the selection criteria and does not reflect 
an astrophysical process. 

North (1993) speculated whether a star can be born 
as a He-weak and "die" as a Si star, or, similarly, whether 
a Si star can become a Cr or a SrCrEu while evolv- 
ing. From Fig. [3] it is obvious that there is a surplus of 
SrCrEu objects near the TAMS. According to the mod- 
els by Schaller et al. (1992), an object with 2.5 Mq, has 
at the ZAMS an effective temperature of about 11150K, 
after 585 Myr it reaches the TAMS with an effective tem- 
perature of 8160 K. During the life time on the ZAMS, the 
stellar atmosphere changes due to rotational and magnetic 
effects as well as due to the diffusion mechanism. These 
effects seem to manifest in the HRD. There is a dramatic 
decline of evolved SrCrEu objects for masses lower than 
2.5 Mq whereas it is the opposite case for masses higher 
than that. There are hardly any evolved Si stars but many 
evolved SrCrEu objects. We therefore conclude that there 
are three different effects: 

1. Si stars change their peculiarity as they evolve to- 
wards temperatures cooler than 10000 K on the main 
sequence 

2. There is an obvious transition between Si, SiCr and 
SrCrEu objects at 10000 K 



3. During the evolution of CP2 stars with masses lower 
than 2.25 Mq, there is a significant decrease of the 
incidence for peculiar objects at 8000 K. 

5. Conclusions 

We have investigated a sample of 182 CP2 stars for which 
homogeneous classifications and accurate parallax mea- 
surements are available. The sample was further divided 
into Si, SiCr and SrCrEu objects, the three main groups. 
For all program stars effective temperatures, luminosities 
and thus ages as well as masses were calibrated with the 
help of standard evolutionary models. For this purpose 
we used data and calibrations within the Geneva 7-color 
and Stromgren uvby/3 photometric system. Special care 
was taken to include the appropriate corrections for the 
spectral and thus flux peculiarities. 

The evolutionary status of CP2 stars is especially inter- 
esting because it still not clear whether these objects show 
their peculiar nature soon after arriving at the ZAMS or 
after about 30 % of their life time on the main sequence. 
Answering this question from the observational side has a 
direct impact on the theories about the origin of the strong 
stellar magnetic fields. One model describes the survival 
of frozen-in fossil fields originating from the medium out 
of which the stars were formed, the other one follows the 
idea that a dynamo mechanism is acting in the interior of 
these stars. 

From the location of our program stars we are able to 
conclude that the CP2 phenomenon occurs already at the 
ZAMS. At least 16% of the investigated CP2 stars have 
relative ages below 20%. We find no significant variation of 
the magnetic field strengths over the whole main sequence. 
The effective temperature and the magnetic field strength 
are the main parameters that determine the kind of the 
chemical peculiarity, i.e silicon - chromium - strontium 
and europium (from hotter to cooler temperatures). 

The difference in the samples of CP2 stars with masses 
above and below 3Mq reported by Hubrig et al. (2000) 
seems to be a selection and thus bias effect of the peculiar 
and normal sample. 

We find a strong decline of evolved SrCrEu objects for 
masses lower than 2.5 Mq, whereas there are few evolved 
Si stars but many evolved SrCrEu objects above this stel- 
lar mass. This can be interpreted as Si stars changing their 
peculiarity as they evolve at temperatures cooler than 
10000 K which results in a transition between Si, SiCr and 
SrCrEu objects at that point. Furthermore, the incidence 
of CP2 stars with masses lower than 2.25 Mq decrease to 
almost zero below effective temperatures of 8000 K. 

These strict observational constraints have to be in- 
cluded and explained by theories that model the stellar 
formation and evolution of stars between 1.5 and 4.5 Mq 
in the presence of strong magnetic fields. 
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